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ABSTRACT. This paper analyze an MIXI/G/1 queue with Bernoulli feed back 
random server breakdowns, Bernoulli schedule server vacation and restricted 
admissibility policy. Customers arrive in batches following compound Poisson 
process and are served one by one according to first come first served basis. The 
service time follows general (arbitrary) distribution. The customer feedback to 
the tail of original queue for repeating the service until the service be successful 
with probability p or the customer departs the system if service be successful 
with probability 1-p. After completion of a service the server may go for a 
vacation with probability 0 or continue staying in the system, to serve a next 
customer, if any with probability 1—0. Unlike the usual batch arrival queueing 
model, there is a restriction over the admissibility of batch arrivals in which 
not all the arriving batches are allowed to join the system at all times. This 
restriction admissibility policy is different for different states of the system. 
The system may breakdown at random in accordance with Poisson process, 
while the repair time follows exponential distribution. We obtain the time 
dependent probability generating function in terms of their Laplace transforms 
and the corresponding steady state results explicitly. Also we derive the system 
performance measures like average number of customers in the queue and the 
average waiting time in closed form. 


1.Introduction 


The research on queueing theory has been extensively developed due to 
a lot of significance in the problems relating with decision making process and it has 
made a tremendous impact in industry and logistics sector apart from its significant 
applications in many other areas like air traffic, bio-sciences, population studies, 
health sectors, manufacturing and production sections etc. Numerous authors have 
put their efforts in studying queueing models by importing many aspects on them. 
Also the analysis of queueing models with Bernoulli feedback had a prominent role 
in queueing theory. Takacs (1963) introduced the concept of feedback in queue- 
ing model. Krishnakumar et al.(2001) studied M/G/1/1 model with feedback by 
considering both regular and optional service. Choudhury and Paul (2005) in- 
spected the M/G/1 system with two phases of heterogeneous service and Bernoulli 
feedback. Choi et al.(2003) analyzed two phase queueing system with vacation 
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and Bernoulli feedback. Badamchi Zadeh and Shahkar (2008) analyzed two phas- 
es queueing system with Bernoulli feedback and Bernoulli Schedule server vacation. 


Most recently, research studies on queues with server breakdown have 
been attracted, as an important area of queueing theory and have been studied 
extensively and successfully due to their various applications in production, com- 
munication systems. Mostly in the queueing literature, the server may be assumed 
as an reliable one, so that it never fails. But in real scenario, mostly servers are 
unreliable, as we often encounter the cases where service stations do not last for 
ever and which can encounter break downs. Such phenomena always occur in the 
area of computer communication networks and flexible manufacturing system etc, 
since the performance of system may be severely affected by server breakdowns. 
Regarding restoration of server as the server has to resume its service, it needs 
a repair process immediately. such systems with break down and repair are well 
worth to investigate from the queueing theory point of view as well as reliability 
point of view. Recently, there are several contributions in the research study of 
non-Morkovian single server queueing system, in which the server may experience 
with break downs and repairs. Thus queueing systems with repairable server has 
been studied by many authors including Avi-Itzhak and Naor (1963), Graver(1960), 
Takine and Sengupta(1997), Tang(1997) etc. 


Another aspect in queueing theory is the study of server vacations which 
play a significant role in queueing models. During the last three or four decades, 
queueing theorists are interested in studying queueing models with vacations elabo- 
rately, because of their applications and theoretical structures in real life situations 
such as manufacturing and production systems, computer and communication sys- 
tems, service and distribution systems, etc. Queue systems with batch arrival and 
vacation time was analyzed by Baba (1986), he derived the queue size distribution, 
waiting time distribution and busy period distribution. The more detailed survey 
about vacation queues can be found in the survey of Doshi (1986). The most re- 
search works which have been done in recent past by some researchers on vacation 
models including Levy and Yechiali (1975), Keilson and servi (1987), Choudhary 
(2002), and Madan(2009) etc. 


Control of queues is also one of the most remarkable and interesting 
areas of queueing theory. In this context, there are many researchers including Rue 
and Rosenberg (1981), Neuts (1984), Stidham (1985), Lee and Srinivasan (1989) 
and Huang and Mc-Donald (1998) dealt with control policies on arrivals in queues 
and queueing networks. One of the queue control is, control on admission of ar- 
riving customers which is an interesting control problem for queues. In some of 
queueing models, queue lengths are controlled by the rejection of some of the in- 
coming arriving customers. Such a control policy is called restriction admissibility 
policy and which has been studied by many authors. Madan and Abu-Dayyeh 
(2002, 2003) proposed restricted admissibility policy on M*¥ /G/1 queueing system 
with Bernoulli server vacation. Madan and Choudhury (2004), Alnowibet and Tad- 
j(2007) and Choudhury (2008) studied many queueing models in this perspective. 
Badamchi Zadeh (2012) studied a batch arrival queueing system with two phases of 
heterogenous service with optional second service and restricted admissibility with 
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single vacation policy. 


Transient state indices are very essential to track down the functioning 
of the system at any instant of time. Thus in this paper, we propose a study on 
transient behavior of a queueing model with feedback subjected to random break 
downs and Bernoulli vacation model under restriction policy for a single server 
batch arrival in put. If a customer feels that, his service is being unsatisfactory 
he require the service repeatedly with probability p or the customer departs the 
system if service be successful with probability 1-p. once the system gets break 
down, it enters a repair process and the customer whose service is interrupted goes 
back to the head of the queue. Also the server may go on vacation according to 
Bernoulli scheduled server vacation, that is, after a service completion, the server 
may go for a vacation with probability 6(0 < @ < 1) or may continue to serve the 
next customer, if any, with probability 1 — 6. The service time and the vacation 
time are generally distributed, while the repair time is exponentially distributed. 
The customers arrive in batches to the system and served one by one on a ” first 
come - first served” basis. Also in our model restricted Admissibility policy has 
been applied that, not all arriving batches are allowed to join the system at all 
times. 


The rest of the paper has been organized as follows: in section 2, the 
mathematical description of our model has been found, in section 3, the transient 
solution of the system has been derived, in section 4, the steady state analysis has 
been discussed. 


2. Mathematical description of the queueing model 


To describe the required queueing model, we assume the following. 


e Let Acidt; i = 1,2,3...be the first order probability of arrival of ’i’? customer- 
s in batches in the system during a short period of time(t,t+dt)where 0 < ci < 
1,5 = 1, > 0 is the mean arrival rate of batches. Unlike the usual 
batch arrival queueing systems, considering the policy of restricted admissibility 
of batches in which not all batches are allowed to join the system at all times. Let 
rı(0 < rı < 1) be the probability that an arriving batch will be allowed to join the 
system during the period of server’s non-vacation period which includes system’s 
service time and repair time also let r2(0 < rg < 1) be the probability that an 
arriving batch will be allowed to join the system during server’s vacation period. 


e After completion of service, if customer received unsuccessful service for certain 
reasons, the customer may immediately join the tail of the original queue as a feed- 
back customer for receiving another regular service with probability p(0 < p < 1). 
Otherwise the customer may depart forever from the system with probability q ( 
= 1 - p). The service discipline for feedback and newly customers are first come 
first served discipline. Also service time for a feedback customer is independent of 
its previous service times. 


e There is a single server which provides service following a general(arbitrary) 
distribution with distribution function B(v)and density function b(v). Letu(x)dx 


4 DR. G. AYYAPPAN AND S. SHYAMALA 


be the conditional probability density function of service completion during the 
interval (x,x+dx] given that the elapsed service time is x, so that 


(1) ulz) = ay 


and therefore 


(x) 


v 
a I ù 
(2) bv) = (ve Jo” “de 
e As soon as a service is completed, the server may take a vacation of random 


length with probability 6 (or) he may stay in the system providing service to a 
customer if any, with probability 1 — 6, where 0 < 0 < 1. 


e The vacation time of the server follows a general(arbitrary)with distribution 
function V(s)and the density function v(s). Let v(x)dx be the conditional prob- 
ability of a completion of a vacation during the interval (x,x+dx] given that the 
elapsed vacation time is x so that 


(3) v(a) = 5 


and therefore 


- E 
(4) v(s) = v(s)e dx 


e The system may breakdown at random and the breakdowns are assumed to 
occur according to a Poisson stream with mean breakdown rate a > 0. Further we 
assume that once the system breakdown, the customer whose service is interrupted 
comes back to the head of queue. 


e Once the system breaks down it enters a repair process immediately. The 
repair times are exponentially distributed with mean repair rate 8 > 0. 


e Various stochastic processes involved in the queueing system are assumed to 
be independent of each other. 


3.Definitions and Equations governing the system 

We let, 

(i)Pa (x,t) = Probability that at time ’t’ the server is active providing service 
and there are ’n’ (n > 0) customers in the queue excluding the one being served 


and the elapsed service time for this customer is x. 
Consequently p,,(t) denotes the probability that at time ’t’ there are ’n’ customers 
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in the queue excluding the one customer in service irrespective of the value of x. 


(ii) Vn (x,t) = Probability that at time ’t’, the server is on vacation with elapsed 
vacation time x, and there are ’n’ (n > 0) customers waiting in the queue for service. 
Consequently V,,(t) denotes the probability that at time ’t’ there are ’n’ customers 
in the queue and the server is on vacation irrespective of the value of x. 


(iii)Rn (t) = Probability that at time t, the server is inactive due to breakdown 
and the system is under repair while there are’n’ (n > 0) customers in the queue. 


(iv)Q(t)= Probability that at time ’t’ there are no customers in the system and 
the server is idle but available in the system . 


3.Transient solution of the queueing model 


The model is then, governed by he following set of differential-difference equa- 
tions. 


(5) 
a a < 
gre taz Pale, t)+(A+u(x)+a)Pa(x, t) = AL — r1)Pa (x, t)+àrı Me C:Pn—i(a,t);n > 1 


i=1 


(6) E + P Bili + (A+ u(x) +a)Polx, t) = A(1 — rı)Po(x, t) 


ot Ox 
(7) 
a a = 
g E t)t+ 5 Vals, t)+(A+v(2))V, (a, t) = AC — re) Vp (£, t)+Are >, CiVn—i(z,t);n > 1 
a) a 
(8) gaVe t) + By Vole: t) + (A+v(x))Vo(z,t) = A(1 — r2)Vo(z, t) 
9) p 
L Rall) = — (A + B)R,(t)+A(1 — r1 )Rn (x, t)+àrı 2 GRn_-i(2, ita | Pr-1(a, t)dx 
(10) © Rolt) = (A+ B)Rolt) + AL — 1) Rol, t) 


dt 
£Q(t) = —AriQ(t) + Ari (1 — ra) Q(t) + r28 Ro(t) 


(11) +rı a Vo(a, t)v(x)dz + r2(1 — oaf P(x, t)u(x)dx 


The above equations are to be solved subject to the following boundary condi- 
tions 
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PON =r- fp Pr(etyu(ade +a f Parle duode} 
0 0 
(12) 
+n f Va+ı(x, t)y(a)da + reBRr+i(t) + Arirecn41Q(t);n > 0 
0 


(13) rıVn(0,t) = r0 | P,(x,t)ulx)dz;n > 0 
0 
Assuming there are no customers in the system initially so that the server is idle. 


(14) Vo(0) = 0; Va (0) = 0; Q(0) = 1; Pa (0) = 0,n = 0,1, 2,... 


Generating functions of the queue length 
The time dependent solution. 


We define the probability generating functions 
P,(2, z,t) = eae 2" P, (x,t) 
P,(z,t) = Dan 2” Pa(t) 
Va(a, 2, t) = ono 2” Vn (z, t) 
Valz, t) = po 2” Vn (t) 


Ra(z,t) = Era 2” Ra (t) 


(15) C(z) = X enz” 


which are convergent inside the circle given by |z| < 1 and define the Laplace 
transform of a function f(t) as 


(16) Hs) = [pleat 


Taking Laplace transforms of equations (5) to (13) 


(17) 


2 Bls) +(s+At+ p(x) +a)Py(z,8) = AC = ri) Palæ, 8) + Ary) Paila 8) 








(18) Ž Ryle, s) h(s +A + p(x) +a)P(z, s) = A0 — rı) P (z, s) 
(19) Žv, s) + (s +A +v(£)) Valz, 8) = A(1 — r2)Vn (z, 8) +Ar2 3 ciVn—ilx, 8) 
(20) OTi, s) + (s +à + v(z))Vo(z, s) = A(1 — r2) (x, 8) 


Ox 
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n—-1 


(21) (s+A+8)Rn(s) = A(1 — 71) Rn(z, 8) +Ar1 2 ciil) ta | Py-i(a, s)dx 


(22) (s +à + B)Ro(s) = A(1 — 11) R, (2, 8) 
(s + Ar1)Q(s) = 1 + Ari(1 — r2)Q(s) + BreRo(s) 


(23) 
+ro(1— Da f Po(a, s)u(a)da +n f Vo(a, s)v(x)dx 


for boundary conditions, 


KEG RS O is fv v AE A f BN s)u(z)ae} 


(24) 





+n f Vas 2, s)v(x)dr+re8 Rn+i(s)+Arir2ACn41Q(s);n >0 
0 


(25) r1V, (0,8) = r0 | P, (x, s)ulx)dz;n > 0 
0 


multiply equation (17) by z” and add (18) implies and using the probability 
generating function defined above. 


(26) 2 P,(0,2,5) H (s+ Ari(1— C(z)) + u(x) 4 a)P, (x, z, 8) =0 


performing similar operations to equations (19) to (22). 











(27) 2 Vale, 218) + (s + Are(1 — C(z)) + v(x))Vq(z, z, 8) = 0 
(28) (s+ Ari (1 — C(z)) + 8)R,(z, 8) =a f P (£, z,8)dx 


For the boundary conditions, we multiply equation(24) by z”+!, sum over n from 
0 to œ and use generating function defined above, we get 


2reP, (0,23) = ra(1=0)(pe +0) | P,(æ,2,s)ule)de+ri f V,(x, z, 8)v(x)dx 
(29) +ro8Rq(z,s) + (1 — sQ(s)) + Arire(C(z) — 1)Q(s) 


Similarly multiply equation (25) by z”and sum over n from 0 to oo and use 
generating function defined above 


(30) riV,(0, z, s) = r0 | P (z, z, s)u(x)dz;n > 0 
0 
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Integrating equation(26) from 0 to x yields 





a = (s+Ari(1—C(z))+a)a J u(t)dt 
(31) P(x, 2,8) = P,(0,z, s)e 0 


where P,(0, 2,8) is given by equation(29). Again integrating equation (31) by 
parts with respect to x yields 














= ee a e ”) 
where 
(33) B(s + Ari (1 — C(z)) +a) = e7 (sari (1-C(2))+0)q B(x) 


is Laplace - Stieltjes transform of the service time B(x). Now multiplying both 
sides of equation (31) by u(x) and integrating over x, we get 


(34) i P (a, z,8)u(x)dx = P,(0, z, s)B(s + Ari(1 — C(z)) +a) 
Using equation (34), equation (30) becomes 


(35) r1V,(0, z, s) = T20P,(0, z, s)B(s + Ari(1 — C(z)) + a) 


Similarly integrate equation (27) from 0 to x, we get 





7 7 (stAre(1—C(z)))a J v(t)dt 
(36) V(x, z, s) = V0, z, s)e 0 


substituting the value of V,(0, z, s) from equation(35) in equation (36) we get 
(37) 
-tantc | v(t)dt 
0 





V,(z, 2,8) = (2) OP, (0, z, s)B(s + Ari(1 — C(z)) + a)e 


Again integrating equation (37) by parts with respect to x 





(38) 
Z (z.s) = (2 \0B 0.2.9)Bls + Ar aita 1— V(s + Are(1 — C(z))) 
aa (2) ero EAE | (s + àr2(1 — C(z)) | 
where 
(39) V(s+r-acly) = f eT (starz -CTAN (x) 
0 


is Laplace - Stieltjes transform of the vacation time V(x). Now multiplying both 
sides of equation(37) by v(x) and integrating over x, we get 
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(40) 
| V,(2, z, 8)v(x)dz = (2) OP, (0, z, s)B(s + Ari(1 — C(z)) + a)V(s+Are(1 — C(z))) 


rı 





Using equation (32), equation (28) becomes 


azP,(0, z, s)[1 — B(s + Ari(1 — C(z)) + a)] 
[s+ Ari (1 — C(z)) + Alls + Ari(1 — C(z)) + a] 


Now using (34), (40) and (41) in equation (29) and solving for P,(0, z, s) we get 


(41) R,(z, 8) = 








fil) fal2)| (4) (1 - 8Q(s)) + àri (Cl) — NQ(s) 
Dr 





(42) P,(0, z,s) = 


where 


Dr = fi(z)folz){z — (1 — 0) (pz + g) BU fi (z)] — 0V (s + Ara — CEMB 
—aßz{1 — Blfr(z)]} 
fi(z) =s4+Ari(1—C(z)) +a 











fo(z) = s +Arı(1-— C(2)) +B 





substituting the value of P,(0,z,s) from equation (42)in to equations (32), (38) 
and (41) 


BO- BEME) (1 = 5Q(s)) + àr (Ci) - NQ(s)| 
Dr 








- fi(2) fale) (2) BAE) a- s8) Aree) - NAs)] ER 





Dr 


az[1 = Bifi(2)]]| (4) (1 = sQl) + Ari (C) - NA(s)| 
Dr 


(45) — -Ra(z,8) = 





where Dr is given above. 


4.The steady state analysis 


In this section we shall derive the steady state probability distribution for our 
queueing model. To define the steady state probabilities, suppress the argument’t’ 
where ever it appears in the time dependent analysis. 

By using well known Tauberian property, 


(46) Liat (s) = sth) 
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multiplying both sides of equation (43), (44), (45) and applying equation(46) and 
simplifying, we get 


_ A(z) -= BlA@))Pri(C((2)) = DA] 











(47) P,(z) A 
(48) V,(2) = PDE) PIA VIN rat -oe - 119 
(49) R,(z) = az|1 — PAA (CUA) = 1)Q] 


where Dr and fı(z) and f2(z) are given in the previous section. Let W,(z) denotes 
the probability generating function of queue size irrespective of the state of the 
system. Then adding (47), (48) and (49) we get 


W,(z) = POC- Bifi Pra (C(@)=-YQ) 








i PA BIAT Ara-C) 


a2(1 - BIA OAr (C) - DQ] 
Dr 
In order to obtain Q, using the normalization condition, as follows 





(51) m 


(52) W,(1)+Q=1 


We see that at z=1, W,(z) is indeterminate of the form 0/0. We apply L’Hospital 
rule in equation (51) 


_ AN QED) (a + 8)[L — B(a)] + Ar2008B(@)E(DE(V) 
dr 
where dr = a8(q + p0)B(a)—An E(D)[(a + B)(1 — B(a))|—0Ar2a6B(a)E(1)E(V) 





(53) W,(1) 





and B(0) = 1,V(0) = 1,-V’(0) = E[V] the mean vacation time. Using equation 
(53) in equation (52) 


_,_ AEU) 1. 1 1 1), 
(54) Q=1 aon” aw + ew 3 Z| + erezinewy} 


and the the utilization factor p of the system is given by 








AET) 1 1 1 1 
E a A 
where p < 1 is the stability condition under which the steady state exists, equa- 
tion(54) gives the probability that the server is idle. Substitute Q from equation 
(54) in equation (50), W,(z) have been completely and explicitly determined which 
is the the probability generating function of the queue size. 
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The average queue size and average waiting time can be obtained as follows 
Let L,denote the mean number of customers in the queue under the steady state, 
then Lg = £W,(z) |z=1,since this formula gives 0/0 form, then we write 


W(2)= a where N(z) and D(z) are the numerator and denominator of the right 


hand side of equation (49) respectively, then we use 





D'(1)N"(1) — N’(1)D"(1) 
2[D’ (1)? 





(56) L= 


where primes and double primes in equation (56) denote first and second derivation 
at z=1 respectively. Carrying out the derivatives at z=1, we have 


(57) N'(1) = AE(1)Q[ri(a + B) — B(a)(Or2aBB(V) — ri(a + B))] 


N"(1) = 2Q—DE(DP {G2 — r?) + Bla? - 3B - O(a + BEV) 





+1008r3B(V?) + B'(a)(r?(a + 8) — 6rireaBE(V)) } 





(58) +AQE(I(I — 1)) {ri(a + B) + B(a)(6r2aBE(V) — ri(a+ 8))} 
(59) 
D'(1) = —AE(D)ri(a + b) + B(a) {ri0B(q + p9) + AE(D)ri(a + B) — OaBr2E(V)} 





D"(1) = 2AE]? {(r? — G+) + Bla) 


[-(4 + p9)ri + Orira(a + B)E(V) — 3080r3E(V?)] 





+B" (a)[-(q + pO)r?(a + 8) — HER + BaBrireE(V)]} 


(60) +AE(I(I — 1)) {=ri(at 8) + B(a)(ri(a + B) — OrgaBE(V))} 


where E(V7) is the second moment of the vacation time and Q has been found 
in equation (54). Then if we substitute the values of N’(1), N” (1), D’(1) and D’(1) 
from equations (57), (58), (59) and (60) in to (56) equation we obtain L, in a closed 
form. 
Mean waiting time of a customer could be found, as follows 


L 
(61) Wi = z 
by using Little’s formula. 
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